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Abstract A precision measurment of inclusive electron scattering cross sections was carried out at Jefferson 
Lab in the quasi-elastic region for 4 He, 12 C, 56 Fe and 208 Pb targets. Longitudinal (Rl) and transverse (Rt) 
response functions of nucleon were extracted in the momentum transfer range 0.55 GeV/c<|g|<1.0 GeV/c. 
To achieve the above goal, a Nal (Tl) calorimeter was used to distinguish good electrons from background 
including pions and low energy electrons rescattered from walls of the spectrometer magnets. Due to a large 
set of kinematics and changes in HV settings, a number of calibrations were performed for the Nal (Tl) detector. 
Corrections for a few blocks of Nal (Tl) with bad or no signal were applied. The resolution of Nal (Tl) detector 
after calibration reached w 3% at E=l GeV. The performance of Nal (Tl) detector was compared with a 
simulation. 
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1 Introduction 

The study of nucleon properties in a nuclear 
medium is an essential objective in nuclear physics. 
The Coulomb Sum Rule(CSR) provides one of the 
cleanest means to study nuclear medium effects on 
the charge response of the nucleons[l|. The Coulomb 
sum S L (q) is given by 



(1) 



G E {Q 2 ) 2 

with Z the atomic number of the nucleus, Q 2 the four 
momentum transfer squared, q the three momentum 
transfer and w the energy loss. After factoring out an 
effective nucleon charge form factor G E {Q 2 ) 2 which 
is an appropriate sum of neutron and proton charge 
form factors and the longitudinal response R L (q,ui) is 
integrated from uj + to infinity where uj + is selected to 
exclude the elastic peak, the S L (q) should approach 
1 as q — > oo for a system of non-relativistic nucleon. 



In this limit S L (q) simply measures the total charge 
divided by Z. In the Fermi gas model the asymptotic 
limit of S L (q) is reached for q > 2k F ~ 500 MeV/c 
where correlations due to the Pauli Blocking effect 
vanish. Since the ratio of R L to R T is small at large 
q, R L has a large sensitivity to the uncertainties of 
the cross sections. To make a precision measurement 
of R L , the uncertainties of cross sections are needed 
to be at 1% level. Therefore, R L was much harder to 
determine with good precision than R T . 

Jefferson Lab CSR experiment (E05-110)Q mea- 
sured the cross sections of quasi-elastic electron scat- 
tering on four different targets ( 4 He, 12 C, 56 Fe and 
208 Pb) at four different scattering angles (15°, 60°, 
90°, 120°) with beam energies from 0.4 GeV to 4.0 
GeV. The standard Hall-A detector configuration in- 
cludes two high resolution spectrometers (HRS)(Hj. 
Each HRS has a Q1Q2-DQ3 magnet configuration 
where Ql, Q2 and Q3 are quadrupole magnets and 
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D is a dipole magnet. For the CSR experiment, both 
HRSs were configured for electron detection. The Nal 
(Tl) detector was installed in the left HRS. We will 
focus on the left HRS. The left HRS detector pack- 
age consisted of two Vertical Drift Chambers (VDCs), 
a pair of plastic scintillator planes, a gas Cerenkov 
counter and a Nal (Tl) calorimeter. The VDCs were 
used to determine particle trajectory. The scintilla- 
tors made the trigger. The gas Cerenkov counter and 
the Nal (Tl) calorimetor formed the particle identifi- 
cation (PID) system. 

2 Nal (Tl) Detector 

This Nal (Tl) calorimeter was first used at Los 
Alamos National Laboratory^ and Brookhaven Na- 
tional Laboratory. The detector was transferred to 
Jefferson Lab for this experiment. The Nal (Tl) 
calorimeter was refurbished and reconfigured into 
three boxes with each box consisting of 90 (10x9 ar- 
ray) blocks. The length, width and height of each 
individual block are 30.5cm, 6.35cm and 6.35cm, re- 
spectively. Because the total length of 30.5 cm was 
11.5 radiation length, an electron with less than 0.55 
GeV could deposit most of its energy in the calorime- 
ter. An electron with energy higher than 0.55 GeV 
would have some energy leaked. Since a few blocks 
have bad or no signal during the experiment, the 
missing energy corrections for the bad blocks were 
studied and corrected in the calibration. The follow- 
ing section will focus on the middle box of the Nal 
(Tl) detector. 

3 Nal (Tl) Calibration 

An electromagnetic cascade generates a shower of 
low energy photons and electron-positron pairs when 
a high energy electron hits the Nal (Tl) calorimeter. 
As the cascade propagates, a large part of the original 
particle's energy is converted to lights, which usuall 
covers several blocks. The light in each block is col- 
lected in a photomultiplier tube (PMT). The output 
signal from the PMT is then digitized with an analog- 
to-digital-converter (ADC). There is a conversion be- 
tween raw ADC values and the total energy deposi- 
tion in the Nal (Tl) calorimeter. To accomplish this 
conversion, one needs to determine calibration coeffi- 
cients for each of the calorimeter blocks. In general 
the electromagnetic cascade is spread over several ad- 
jacent blocks, the output signal must be integrated 
over the entire calorimeter volume to obtain the total 
detectable energy. If the calorimeter has been cali- 



brated properly, the total deposited energy E should 
be proportional to the incident particle's energy (or 
momentum p). 

3.1 Calibration event selection 

To obtain good calibration coefficients, an elec- 
tron sample needed to be selected. This was accom- 
plished by selecting a run with less background from 
pions. Runs in the quasi-elastic electron scattering- 
settings were used to do the calibration. The follow- 
ing tight cuts were applied to select electron sample: 
(1) An event reconstruction in the spectrometer de- 
tector package was successful; (2) Only one track 
had been reconstructed by the VDC system; (3) The 
event was identified as an electron with a tight Gas 
Cerenkov cut; (4) The event was in a good accep- 
tance region of spectrometer; (5) The event was in 
the central region of the block being calibrated. 

3.2 The method for determination of Nal 
(Tl) calibration coefficients 

The calibration coefficients are defined to trans- 
form the ADC amplitude of each block into the en- 
ergy deposition of the electron in this block. Since the 
Moliere radius of Nal (Tl) is 4.8 cm[5], the incident 
particle's energy is deposited in the 9 blocks when it 
hits the central block (i.e. blk5 in Fig. 1). The basic 
calibration cell is set to be 9 blocks. A linear mini- 
mization method is used to determine the calibration 
coefficients. The Chi-square minimization function is 
defined as follows: 

N / 9 \ 2 

X t = y E[EL- y EO k A{\ , (2) 

i=o \ k=o / 

where j is the index of the selected calibration events 
and k is the index of the Nal (Tl) blocks. A{ is the 
amplitude in the k-th Nal (Tl) block for the j-th 
event. E kin is the scattering electron energy 5 is 
the calibration constant for the A:-th block. 

^ = ^t(EL-±C k Ai\ , (3) 

where i varies between and 90. \ 2 IS minimized 
when the above quantity is set to zero. It leads to: 

E ( c * (l>i4 J J =f>L4, ( 4 ) 

k=0 \ \j=0 ) J i=0 

The linear equation can be summarized in matrix 
form ... 

MC = E, (5) 
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where C and E are denned as vectors 



/ Co \ 



C = 



(6) 



E = 



(7) 



and the matrix elements are given by: 

JV 



(8) 



At the end, the calibration coefficients are ob- 
tained by inverting eq. (5): 

C = M' 1 E, (9) 

3.3 Missing energy correction for bad blocks 

To find an average value to correct for the missed 
energy in the neighboring block which was bad, 9 
good blocks were selected for studying the amplitude 
ratio of adjacent block to the central one. After ob- 
taining the relationship, the missing energy of bad 
blocks could be corrected back for calibration. A few 
small circle cuts in the region of central block were set 
to obtain the amplitude ratio for scattering electron 
with p= 539 MeV/c at 60°. 
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The amplitude ratios of adjacent blocks to the cen- 
tral one (|x| < 0.03m, \y\ < 0.03m) were fitted to 2nd- 
order polynomials in two dimensions. The following 
relations were obtained: 

i?l = 0.01473-0.09517x + 7.359x 2 + 0.1848?; + 7.51?/ 2 , 

(10) 

R2 = 0.01402-0.3942a;+8.537a; 2 -0.0014132/-2.994y 2 , 

(11) 

i?3 = 0.01399 - 0.1175x + 4.114x 2 - 0.1725?; + 4.55jy 2 , 

(12) 

i?4 = 0.02414-0.08904x+21.29:r 2 +0.8816y+23.54y 2 , 

(13) 

R5 = l, (14) 

R6 = 0.02423-0.001044x+19.357a; 2 -0.5776y+9.302y 2 , 

(15) 

R7 = 0.008895-0.09234a;+6.1917a; 2 -K).09544j/44.266 2 / 2 , 

(16) 

R8 = 0.01739+0. 7615x+22. 78x 2 +0.01578y+5.965y 2 , 

(17) 

R9 = 0.01359 + 0.1170x + 7.258x 2 - 0.2098?/ + 7.97y 2 , 

(18) 

where Rl, R2, R3, R4, R5, R6, R7, R8 and R9 are the 
amplitude ratios of each adjacent block to the central 
block, respectively. The x and y represent the vertical 
and horizontal directions shown in Fig. 1 Since the 
amplitude ratios of adjacent blocks has a momentum 
dependence, different math forms were taken for dif- 
ferent momentum settings for which calibrations were 
performed. 

4 Checking the calibration results 

The energy deposition E of an incident electron 
in the calorimeter detector could be calculated with 
the calibration constants C by the following formula: 



E-- 



A, 



(19) 



Fig. 1 9 block scheme for NaT (Tl) calibration 



where i was the number of Nal (Tl) detector block. 
The E/p of electrons should be around 1 after the 
calibration correction where E was the total deposit 
energy calculated from Eq. (19) and p was the mo- 
mentum. The E/p plot before calibration is shown in 
Fig. 2 at 120 MeV of scattering electron. The plot 
after calibration correction is shown in Fig. 3. 
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Fig. 2 The E/p plot before calibration for a 
non-calibration run 

The E/p plot for electrons should be a Gaussian dis- 
tribution plus a tail. The width of the Gaussian dis- 
tribution represents the detector's resolution. Fig. 3 
shows a reasonable spectrum after calibration correc- 
tion. 
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Fig. 3 E/p plot after calibration for the electron 
peak of a non-calibration run 

The events at the low momentum tail were electrons 
with energy leakage, electrons that scattered the walls 
of the spectrometer and secondary productions, in- 
cluding electrons and hadrons. The energy resolution 
of E/P in Fig. 3 is about 9.7% for 120 MeV at this 
HV setting and the best resolution we can reach is 3% 
for 1 GeV. Due to a large set of kinematics and the 
changes of the HV for Nal (Tl) detector during the 
data-taking period, different sets of calibration con- 
stants were needed for this experiment. A total of 40 
sets of constants for production runs were obtained 
for this experiment. 

5 Nal (Tl) simulation for background 
analysis 

A simulation using SNAKE Q and GEANT30 
was performed when the experiment was proposed[2j. 
We used that to generate an input electron sample 



for the Nal (Tl) detector. GEANT4 was used for Nal 
(Tl) detector simulation. 

5.1 SNAKE and GEANT3 simulation 

The background generated by the interaction of 
electrons with the inner walls of the spectrometer 
magnets was studied with a Monte-Carlo simulation. 
The simulation is based on a ray-tracing program, 
SNAKE. In the original version of SNAKE, an elec- 
tron hitting the internal boundaries of the spectrome- 
ter was considered lost. In the modified version of the 
simulation program, the electron was studied further 
with a GEANT3 simulation for one of two possibil- 
ities: (a) scattering off the wall; (b) generation of 
secondary particles from an interaction with the wall 
material. 




Fig. 4 The position and direction of electron 
sample generated by SNAKE and GEANT3 at a 
momentum setting of 120 MeV 

Then the rescattered electron or the secondary parti- 
cles were re-inserted into the SNAKE simulation and 
were traced to the focal plane. Since particles hitting 
the walls of Ql, Q2 could not make to the focal plane 
and the ones hitting the dipole have small probabil- 
ities to reach the focal plane in the SNAKE simula- 
tions. So the simulation was focused on the interac- 
tion of electrons with the walls of the Q3 magnet. Due 
to the proximity of the Q3 magnet to the focal plane, 
electrons bouncing off the surface of the Q3 magnet 
would have a higher probability of survival than those 
bounced off the other magnets. The result of simu- 
lation shows that the background generated in this 
process is about 1.2% of the clean events at a spec- 
trometer momentum setting of 120 MeV/c, the low- 
est momentum setting among our kinematics. A few 
background events with energy comparable to clean 
events came from a single, large angle scattering on 
the surface of the Q3 magnet. With a tight cut on the 
position on the focal plane, about 80% of the back- 
ground events were eliminated. This is in agreement 
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with results from an independent analysis [8]. The re- 
maining background events after the focal plane po- 
sition cut can be eliminated by an independent en- 
ergy measurement such as Nal (Tl) calorimeter. The 
position and momentum direction of electron sample 
before entering Nal (Tl) detector which were gener- 
ated from the SNAKE and GEANT3 simulation were 
shown in the Fig. 4 (a) and (b), respectively. The 
electrons reflected by the Dipole and Q3 contributed 
about 0.24%, 1.2% background, respectively. 




10 3 
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Fig. 5 The energy distribution of electron sample 
generated by SNAKE and GEANT3 at a 
momentum setting of 120 MeV 



The fig. 5 shows the energy distribution for the elec- 
tron sample. The good electron peak (dash line) is at 
120 MeV and the background (solid line) spreads out, 
which are electrons scattered by Dipole and Q3 mag- 
nets. The electrons were re-inserted into GEANT4 to 
simulate their behavior in the Nal (Tl) detector. 



5.2 Nal (Tl) GEANT4 simulation 



The geometry of the Nal (Tl) detector is shown in 
Fig. 6. The parameters used in the GEANT4 simu- 
lation for Nal (Tl) properties were obtained from the 
manufacture [9( . The trajectories of electrons and 
photons are shown in Fig. 6 . Since the input elec- 
tron energy was 120 MeV, the total energy of the 
electron was absorbed by Nal (Tl) calorimeter. 




Fig. 6 Schematic of GEANT4 simulation for Nal 
(Tl) detector at a momentum setting of 120 MeV of 
scattering electron; Red line (electron), Green line 
(photon). 

5.3 A comparison between data and simula- 
tion for 120 MeV and 539 MeV 

Fig. 7 shows, in log scale, the distribution of en- 
ergy deposition by electrons in the Nal (Tl) detector 
for simulation result(dash line) and data(solid line), 
respectively. Since the simulation did not contain 
spectrometer intrinsic resolution, hence the simula- 
tion only reproduces the low energy part of Gaussian 
distribution after electron deposited energy in Nal 
(Tl) middle box. Becasuse the low energy background 
part was the main issue, in our case the asymmetry 
of simulation result has no effect on the background 
analysis. 




100 150 200 

Energy distribution (MeV) 

Fig. 7 Comparison between data (solid line) and 
simulation (dash line) at a momentum setting of 120 
MeV 

The good match between simulation and real data in 
low energy part of distribution was shown in Fig. 7 
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for 120 MeV of scattering electrons. From simulation 
result, the contamination from scattering off walls of 
the Dipole and Q3 was about 0.31% and good elec- 
tron cut efficiency was 99.9% after adding a cut on 
the energy distribution at 50 MeV for this kinematic 
setting. 




Energy distribution (MeV) 

Fig. 8 Comparison between data (solid line) and 
simulation (dash line) for 539 MeV 

The fig. 8 shows a comparison between the simula- 
tion and data for scattering electrons with momenta 
of 539 MeV. With a cut applied at 150 MeV, the 



residual contamination from surface scattering from 
the Dipole and Q3 was about 0.029% and good elec- 
tron cut efficiency was 99.9%. 

6 Conclusion 

The Nal (Tl) detector has been calibrated for the 
Coulomb Sum Rule experiment. Corrections were 
applied to the missing energy due to a few inefficient 
detector blocks. Because of the large set of kinemat- 
ics and the changes in HVs, 40 sets of calibration 
constants were obtained. The energy resolution of 
the Nal (Tl) detector reached ^§ » 3% for 1 GeV 
electrons. We also did simulation to study the back- 
ground due to re-scattering from the inner walls of 
the Dipole and Q3 for the spectrometer momentum 
settings of 120 MeV and 539 MeV. The contamina- 
tion was about 0.3% and 0.03% when cuts at 50 MeV 
and 150 MeV were applied for the momentum set- 
tings of 120 MeV and 539 MeV, respectively (see Fig. 

7 and Fig. 8). With the same cuts, good electron cut 
efficiency was 99.9% for both settings. 
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